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Abstract 



We study the diffusion of classical hard-core particles in disordered lattices within 
the formalism of a quantum spin representation. This analogy enables an exact 
treatment of non- instantaneous correlation functions at finite particle densities in 
terms of single spin excitations in disordered ferromagnetic backgrounds. Appli- 
cations to diluted chains and percolation clusters are discussed. It is found that 
density fluctuations in the former exhibit a stretched exponential decay while an 
anomalous power law asymptotic decay is conjectured for the latter. 



PACS numbers: 02.50.-r, 75.10.Jm, 82.20.Mj, 05.50. +q 



1 Introduction 



Processes involving classical particles constrained to diffuse and interact stochas- 
tically on discrete substrates is one of the oldest schemes for the study of a rich 
variety of relaxation phenomena |2[. Although from a conceptual standpoint such 
processes give only a phenomenological level of understanding, they do provide the 
type of behavior necessary to describe irreversibility in many-body systems which 
otherwise, would be difficult to derive from first principles. Among the main tech- 
niques to obtain the time probability distribution of these processes, possibly the 
master equation (MEQ) approach || is the more directly related to physical con- 
cepts. Despite its apparent simplicity however, it generally gives rise to an infinite 
hierarchy of coupled equations whose solutions become quite involved to elucidate, 
particularly at large times. In a concerted effort to remedy this situation, quantum 
field theory methods have recently regained new impetus in the study of the dynam- 
ics of non-quantum- mechanical-objects @]. Basically, the underlying idea is that the 
MEQ resembles a time-dependent Schrodinger equation in a pure imaginary time 
with probability distributions playing the role of wave functions. Thus, by studying 
the field theory associated to the MEQ a formal solution, in some cases exact, of 
many-body probability distributions can be found explicitly ||. 

As a contribution to the current momentum of this approach, in this work we 
shall revisit the problem of hard core particles diffusing in quenched disordered lat- 
tices emphasizing the formal analogy with its quantum counterpart. The avoidance 
of double occupancy accounts for the essential aspects of both inter-particle inter- 
actions and mobility whereas frozen- in disordered bonds {J r , r '}, modeling random 
hopping rates between locations {r, r'}, is essential to understand the slowing down 
of transport properties in a vast family of inhomogeneous and glassy systems [I], || . 

Following the thread of ideas initiated in Q and developed subsequently by many 
authors [|, |8| here we attempt to further advance the subject forward by means of a 
(pseudo) spin- 1 description in which spin up or down at a given site corresponds to 
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particle or vacancy, say, at that location. The MEQ is then equivalent to the action 
of a quantum spin " Hamiltonian" encompassing the original transition probability 
rates. Interestingly, for the case of symmetric diffusion of hard core particles the 
resulting Hamiltonian reduces to an isotropic Heisenberg ferromagnet whose full ro- 
tational symmetry applies irrespective of the hopping disorder. Although this contin- 
uous symmetry is ultimately imposed by conservation of probability throughout the 
Brownian process, its exploitation is not evident without appealing to the quantum 
spin Hamiltonian analogy 0. Particularly, this formalism becomes advantageous 
in analyzing spontaneous fluctuations of the steady state such as non-instantaneous 
density-density correlation functions (e.g. structure factors and scattering func- 
tions). As we shall see, the use of elementary quantum mechanical considerations 
along with selection rules based on the conservation of total angular momentum 
enable an exact treatment of these many-body correlations in terms of single spin 
excitations. 

The layout of the paper is organized as follows. In Section 2, we recast the 
MEQ of the disordered Brownian process as a physical realization of the Heisenberg 
ferromagnet. Section 2 A and 2 B treat in turn the structure of the steady state and 
its non-instantaneous two-point correlations. Section 3 discusses autocorrelation 
functions in a variety of scenarios. An exact treatment of these functions in diluted 
chains is given in Section 3 A, where a stretched exponential asymptotic dynamics 
exhibiting a scaling regime is found. Section 3 B continues with a discussion of 
diluted lattices in higher dimensions. Exploiting scaling pictures for the frequency 
and the density of states of single-excitations on percolation clusters at criticality 
[ jH| , we are led to suggest a slower or anomalous diffusive kinetics. We end the 
paper with Section 4 which contains our conclusions, along with some remarks on 
extensions of the present work. 
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2 Master equation and spin representation 



A central assumption underlying most phenomenological stochastic models is that 
the actual non-equilibrium dynamics of real systems can be approximated by a 
discrete Markovian process and, therefore described by a MEQ. The latter governs 
entirely the time evolution of the probabilities P(s,t) of finding the system in a 
certain configuration \s) at time t. If W(s — > s') denotes the rate or transition 
probability per unit time at which configuration \s) evolves to \s') , the MEQ reads 

8 t P{s, *) = £[ W(s' - s) P(s', t) - W(s - s') P(s, t) ] . (1) 

s' 

Assuming the basis vectors | s ) form an orthonormal set, it is useful to describe the 
ensemble averaged state vector of the system at time t as | P(t) ) = J2 S -P( s > t) \s). 
Thus, starting from a given probability distribution | P(0) ) it can be readily checked 
that the formal integration of Eq. (JJ) yields a state vector solution of the type 

|P(t)> = e- fft |P(0)>, (2) 

where the matrix elements of the transition operator H (or " Hamiltonian" ) are 
constructed as NIT 



s 



'\H\s) = -W(s^s') , s'/s, (3) 



(s\H\s) = W(s->s'). (4) 

The steady states of our stochastic processes correspond to the ground states 
of H of each subspace within which the dynamics takes place, all with the zero 
eigenvalue. Any eigenvalues with positive real parts A > correspond to decaying 
states with lifetime 1/A. Due to probability conservation, clearly every column of 
H adds up to zero, i.e. J2 S ' ( s '\H\s) = . Thus, in passing and for future reference 
it is worth pointing out that a left steady state can be immediately obtained as 

<^l=£(s|, (5) 



which evidently satisfies (ip \H = 0. Also, notice that (ip \P(t)) = 1 Vt. 

Turning to the analogy between diffusion of hard-core classical particles and 
quantum spin-1/2 systems we now represent a particle or vacancy, respectively, at 
site r by spin up or down, i.e. by m r = +1,-1 where m r is an eigenvalue of the z- 
component (say) of the Pauli operator <?(r) for site r. Clearly a typical configuration 
\s) in a general lattice with N locations can be characterized by a direct product 
| m ri , ... , m rjv ) of all the site spinors. So the set {W(s — > s')} = { J rr '} of transition 
hopping rates between connected sites r, r' (e.g. nearest-neighbors (r, r')), then 
describes the following spin-exchange process 

J , 

I . r r | . 

| m ri ... m r = 1 ... 77V = — 1 ••• ^rjv) < > I "In ••• m r = —1 ••• "V = 1 ... mr N ) ■ 

Here and in the following the set {J rr '} can be arbitrarily substitutional^ dis- 
ordered. To construct the associated " Hamiltonian" H whose terms keep proper 
track of both probability conservation and the spin exchange process, it is useful 
to cast the discussion in terms of spin raising and lowering operators <7+, cr~. The 
off-diagonal part of H exchanging configurations as schematized above, is therefore 
given by 

E E W) (AH]*) (S\ = - E ^r' + h.C.) , (6) 

* s'^s (rr'> 

whereas conservation of probability, i.e. Eq.(^), requires the emergence of a diagonal 
part of the form 

E I s ) ( S \ H \ S ) ( s \ = E J rr> [ n r (1 - n r ') + (1 - n T ) h r , ] , (7) 

s (rr'J 

where n r = a+ a~ denote occupation number operators. These terms basically count 
the total number of ways in which a given configuration \s) can evolve to different 
states \s') through a particle hopping attempt between nearest neighbor pairs (rr'), 
weighting each accessible attempt with rate J Tr >. This yields precisely the right hand 
side of Eq. (0). Thus, in terms of usual spin-1/2 Pauli matrices a = (a x , cr y ,(r z ), by 
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virtue of Eqs. and (]7|) the evolution operator H reduces finally to a Heisenberg 
ferromagnet 

#=-" E Jrr' - 1) . (8) 

(rr'> 

Interestingly, due the isotropic nature of the interactions, the stochastic dynamics 
leaves invariant all the components of the total angular momentum S = | J2r 
namely [ H, S ] =0 irrespective of the disordered background of hopping rates J r r / . 
As a consequence of the full rotational symmetry, the calculation of spontaneous 
density fluctuations are simplified remarkably (Section 2 B). In preparation for the 
analysis of those functions, let us first examine the form of the steady state at finite 
particle densities. 

A. The steady state 

Evidently, the fully jammed ferromagnetic state | F ) = | S , S z ) with total spin 
S = N/2 and total magnetization S z = N/2 is a steady state (SS). Since [H , S~ ] = 
we can generate normalized SS | ip m ) with | S , S z ) = \N/2, N/2 — m), i.e. having 
m- vacancies and particle density p = 1 — m/N, and such that H\ip m ) = by 
applying m-times the lowering operator S~ = J2r = S x — i S y , namely 

\^ m ) = a m (S-) m \F), (9) 

where the normalization factor is a m = ^ (N — m)\/N\ . 

Notice that \ip m ) is an equally weighted linear combination of all permissible 
Q = (^j configurations with (N — m) particles. In particular, 

(ip m \ = -!=(j\, (10) 



which is in line with the more basic observation that, aside from normalization pref- 
actors, left and right SS should coincide since Eq (§) involves a self-adjoint evolution 
operator, namely W(s —>■ s') = W(s' —>■ s). From this latter observation, it is 



worth pointing out that detailed balance arises immediately in Eq. (|l|) given that 

p(s) -> i/n v \ s ). 

The structure of these equilibrium states is rather trivial regardless of the hopping 
disorder. This follows by noting that 

| An ) = E" - E" "7* - ^ l F >> (11) 

Jl< •••<Jm 

where the index ordering denotes sums over all reachable different configurations and 
the double prime restricts the sums to rj 7^ p, q. This yields (^j^ 2 ) non-vanishing 
configurations and therefore 

N-2\ 

m ) A rn\ ( m \ 2 



> m |n p n q |^ m ) = -^- = ^i-^J ^- jyrr) p > ( 12 ) 



so density-density correlators decouple in the large system limit. A similar reasoning 
holds for many-point correlators thus yielding spatially uncorrelated SS. Henceforth, 
we address our approach to the correlations between fluctuations that occur sponta- 
neously at differenttim.es which as is known |]12 |, are closely related to the relaxation 



dynamics governing non-equilibrium regimes. 

B. Non instantaneous steady state correlations 

We will be especially interested in the calculation of non-instantaneous joint 
probability distributions ( A(t) £>(0) ) of quantities A, B such as local densities or 
local density correlations, i.e. diagonal operators in the particle or a z representation. 
Here, the brackets indicate an average over histories up to a to time t starting from a 
given probability distribution |P(0)). More precisely, this can be expressed in terms 
of the following discrete path integration 

(A(t)B(0)) = (At) n - 1 J2-J2 p ( s i,0)(s 1 \B\s 1 )W(s 1 ^s 2 ) 

si s„ 

x W(s 2 -> s 3 ) ... W(s„_i -> s n ) (s n \A\s n ) , (13) 
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where the sums run over all possible states \sj), and AtW(sj^i — > Sj) denotes the 
probability of evolving from to \sj) in a single elementary time step At = 

t/(n — 1) |13[. However, by construction [Eqs. (|3|) and (§)], this latter probability 
is given by ( s^- 1 1 — At if | Sj_i ) . Because ^4, £> are taken diagonal in the particle 
representation, 

P(s,0)(s\B\s) = {s\B\P(0)), 

(s\A\s) = ($\A\s), (14) 



so, recalling that (1 — At H) n — ► e in the limit of a large number of steps [13 
Eq. (ID yields 

(*4(t)£(0)) = $ \Ae- Ht B\P(0)). (15) 

Since [ H, S ] = for any disorder realization, Eq. (|15|) allows for a systematic treat- 
ment of spontaneous density fluctuations. In terms of the local density fields referred 
to above, this requires consideration of non instantaneous particle-particle correla- 
tions C r>r '(t) = (n r (t) 7v(0) ) occurring in the SS distribution | P(0) ) = with 
p — 1 — m/N . As usual, it is convenient to choose a basis {[A) = | S , S z , A )} of 
common eigenstates of S 2 , S z and H with eigenvalues S(S + 1) , S z and A respec- 
tively. Inserting an orthogonalized set of |A) states in Eq. (0), by virtue of Eq. ([TOD 
it follows that in the m-down spin subspace C r y(t) can be written as 

C r y{t)= e- xt (tp m \h r \A) (A|rv| Vm), (16) 

levels A > 

where we have restricted the sum to eigenstates |A) different from | ip m ) (which 
has vanishing eigenvalue A = 0), to subtract the time independent correlations, i.e. 
(n r )(n r /) = p 2 . All eigenvalues A are real and positive definite, as by construction 
if is a self adjoint stochastic operator. 

Selection rules can now be applied to simplify remarkably the calculation. Noting 
that the Pauli matrices a z are tensors of rank one, since n r = (1 + <j z )/2, the 
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Wigner-Eckart theorem ensures non- vanishing matrix elements (ip m \ n r \A) only if 
the total spin S of \ij) m ) and |A) differ by or 1, i. e. S = N/2, N/2 — 1 ; in either 
case with S z = N/2 — m . To identify the relevant |A) states, we first compute 
the total spin S of a single spin excitation \<f\) = ^ T ip\(r) a~ \ F) . Since S 2 = 
S z S z + \ { S + S~ + S~ S + ) , it can readily checked that 

f(f-l)f | pa) if A > , 
S 2 |pa>= ' (17) 

[f (f + l)|p A > if A = , 

for which the wave function should satisfy ZlrPA(r) = y/~N 5\ t o . Therefore, the 
wanted |A) states giving non-zero matrix elements in Eq. ( jl6|) are essentially rotated 
versions of the above single spin excitation states. Specifically, recalling that the 
angular momentum algebra imposes 



S~ \S, S z ) = y/(S + S z ) (S - S z + 1) \S, S z -1) , 



we can generate a normalized SS in the m- down spin sector by applying (m — 1)- 
times the lowering operator S~ to the spin excitation | ip\ ) , namely 

|A)=/3 m (S-) m - 1 |pA), (19) 

where the normalization factor arising from the products generated by each appli- 
cation of S~ turns out to be 



0n 



(N-m-l)\ 

1 < m < N - 1 . (20) 



\\ (N -2)\(m-l)\ 

Thus, for |A) ^ \ip m ) these are all the linearly independent states contributing to 
Eq. ([16]), having total spin S = {N/2 — 1) , total magnetization S z = {N/2 — m) , 
and A > 0. So, we are now left with the calculation of ( ip m \ h T \ A ) . This matrix 
element is expanded as 

(ip m \h r \A) = Pm(ipm\n r {S') m ' 1 \^\) 

= f3 m Y.Vx{p)(^ m \n r {S~) m - 1 a p \F). (21) 



To go a further step in our analysis notice that 



n r (S-) m - 1 a p \F) = (m-l)!(l-<f rjP ) 

x £"£"-£" a- x ...a- m _^\F), (22) 

jl< J2< ...<Jm-l 

where the double prime restricts the sums to vectors rj 7^ r, p . Thus, for p 
there are ( terms contributing equally to ( ip m \ h r | A } . Hence, by taking 

into account the wave function constraint J2 r V 9 a(i") =0 so as to ensure a total spin 
S = N/2-1 (A > 0), we obtain 

(-0 m |n r |A) = (m- l)\m\ ( a m p m ^2 (px{p) {1- S TtP ) 



m (N — m) , . . . 

^¥W^ W ' (23) 

where a m is taken as in Eq. (^). Therefore the particle-particle correlations of 
Eq. (|T6|) are finally given by 

^(t) = ^ (1 _~ P) E e- Ai VA (r)#'), (24) 

single A levels 

which at most requires the evaluation of single spin-wave (Bloch) excitations in 
disordered ferromagnetic backgrounds. This is the main result of this section. Cer- 
tainly, the usefulness of the quantum spin operational formalism developed so far is 
subject to the knowledge of such excitations. However, much is known about their 



properties and density of states |10], [141 , an d in fact this enables us to obtain explicit 
results, particularly for large time asymptotic regimes (Section 3). A similar rea- 
soning for non-instantaneous i^-point correlators would involve |A) states with total 
spin S = N/2 — M. , i.e. states with Ai- interacting magnons (Ai = 1, ... In 
this more general case however, analytical progress seems difficult given the scarcity 



of exact results in d > 1, even for ordered substrates [15 
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3 Autocorrelation functions of diluted lattices 



A subcase of Eq. fl24[) entailing particular interest and which is related to the proba- 
bility of hard core random walkers returning to the origin, is that of autocorrelation 
functions C(t). This corresponds to averaging ( p4[) over all locations r = r'. For 
isotropic structures, it is plausible that the number of sites a random walker (par- 
ticle) has visited after a time t is proportional to the volume R d {t), where R is the 
root mean square displacement in a substrate of dimensionality d (either Euclidean 
or fractal). Thus for densities p <C 1/-R d , the probability of returning to the ori- 
gin should in principle scale as C(t) oc R~ d (t). Departures from normal or Fickian 
diffusion (R oc t 1 ^ 2 ) in disordered lattices can therefore be studied by both, i) eval- 
uating the asymptotic behavior of autocorrelations and , ii) averaging the hopping 
distribution. 

Autocorrelations are in turn closely related to the density of states in the sub- 
strate. In fact, due to the normalization of the single level eigenf unctions, after 
averaging over the sites (origins) of a given sample, in the large iV-limit clearly 
Eq. fl24l) reduces to 

C{t) = j* E (25) 

single A levels 

where A p = p(l — p) . So, the form of the density of states for A — * determines 
entirely the large time asymptotic behavior of the autocorrelation function for hard 
core particles at finite densities. 

Before continuing and in preparation for the analysis of autocorrelations in disor- 
dered scenarios, it is worthwhile to pause and consider the finite size scaling regimes 
emerging asymptotically from Eq. (p5|) in regular ordered situations ({J( r r'}} = J)- 
Let us consider briefly a d -dimensional hypercubic lattice with N = L d sites and 
periodic boundary conditions. It is well known that the corresponding A levels are 
2 J £j=i (1 ~~ cos ^ 1 ) with rij — 0, • • • L — 1. From Eq. ( P3|) it is straightforward 
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to check that the autocorrelation function of such finite system factorizes as 



C(t) = A, 



,-2dJi 



L d 



L-l 



e 2 J t cos 2™ 



ra=0 



(26) 



where the last term takes care of the cancellation of the A = contribution. Thus, 
it can be readily verified that there is a scaling regime L — > oo, Jt — > oo for which 
the autocorrelation scales with a universal function T{t) so long as r = Jt/L 2 is 
held constant, namely 



C(L,t) 
^(r) 



(Jt)- d / 2 ^(r), 
A P r d/2 \4(Q) 



(27) 



n=l 



where q = e~ 4n T and ^3 (q) is a Jacobi theta function of the third kind |16| . Hence, 
the typical size of the system sets the time scales (t oc L 2 ) for which the dynamics 
becomes diffusive. Also, recalling that lim ^fr $ 3 (g) = the autocorrelations 
of the infinite system result in long time diffusive tails (as they should), whose 
asymptotic kinetics is given by 

.4 



C{t) 



(Jt) 



-d/2 



(2? 



( 4vr )d/2 

Turning to more general situations, of interest is that of bond- (hopping-) 

diluted lattices. More specifically, we will focus attention on the following binary 
hopping probability distribution 



P(J')=p8j,j> + (l-p)6j,, Q 



(29) 



defined uniformly throughout all lattice bonds. For short-range hoppings, say e.g. 
nearest-neighbors, a purely geometric effect, known as the percolation transition, 
arises from the configurational aspects caused by dilution. We address the reader 



to Ref. |17| for a detailed introduction to this subject which has been found useful 



to characterize a rich variety of diluted systems including spin systems []18||. Here 
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we just mention that there exists a critical concentration p c below which only finite 
clusters exist and above which a cluster spanning the (infinite) lattice is formed. 

Of course in d = 1 any bond removal disconnects the lattice, so p c = 1 . However, 



as we shall see in Section 3 A, as in the regular case [Eq. (27) ], autocorrelations also 
exhibit an asymptotic scaling regime for finite dilutions but through a dilution- 
dependent scale of length. In higher Euclidean dimensions, extensive research [|K], 
[TJ], [H], [TJ| has led to the conclusion that percolation clusters are statistically self- 
similar at the transition, and so are random fractals. They are characterized by 
various scaling dimensions (fractal, spectral) for the various processes (mass, density 
of states) measurable on the fractal. Thus, exploiting scaling results for the spectral 
dimension, in Section 3B we shall discuss the implications which the resulting density 
of states imposes on the large time average behavior of Eq. (p5|). 

A. Diluted chain 

In the one-dimensional case for < p < 1 the distribution (|2~9"D breaks the N- 
site chains into a collection of finite disconnected segments, each having a number 
1<L<A — lof consecutive bonds, or L + 1 correlated sites. So, the probability 
W(L) to find a chain with L consecutive non- vanishing transition rates and free 
boundaries is independent of the chain location and is given by 

' (1 -p) 2 p L if 1 < L < A -2, 
W(L) = < (1 -p)p L if L = N-1, (30) 

, p L if L = N (periodic chain). 

On the other hand, as shown in Section 2, for a given segment of (L — 1) bonds, it 
is sufficient to consider the corresponding evolution operator Hl within the single 
spin excitation sector. From Eq. (pi), it can be readily checked that Hl can be recast 
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in terms of the following L x L tridiagonal form 



H T =2J 



1 -1 

-1 2 

'■• 

••• 





-1 



-1 





-. 

2 -1 
-1 1 



(31) 



with eigenvalues 



/ fin \ 

\ L {n) = 4 J f 1 - cos — J , n = 0,1, 



1. 



(32) 



Assuming a constant particle density p per chain, the overall autocorrelations must 
be averaged over both histories (time) and samples (hopping disorder). Thus, the au- 
tocorrelations arising from all possible distributions of disconnected segments reads 



N-l L 

(C(«)> = A p £ W(L) £ 

L=l 



(n)t 



(33) 



n=l 



where A p is taken as in Eq. Here the first sum, denoted by the brackets 

on the left hand side, runs over segments of L consecutive bonds and takes into 
account the hopping disorder, whereas the second sum runs over eigenvalue levels 
Ax+i(n) > 0, 1 < n < L, and carries out the time dependent average. We are 
especially interested to elucidate the asymptotic behavior of Eq. ( pB]) for arbitrarily 
large times and finite dilution regimes < p < 1 . To this aim, it is useful to 
introduce the natural length scale emerging in the problem, namely the percolation 
correlation length 



hip 



< p < 1 , 



(34) 



which measures the effective distance between missing hopping rates. On general 
grounds, this scale can be regarded as the mean distance between two sites belonging 
to the same cluster which, in higher dimensions, diverges as £ ~ (p — p c )~ u with a 
universal exponent v depending solely on the space dimensionality ITT" 
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Now, from Eqs. (j30|) and ( p3|) it follows that in the thermodynamic limit iV — > oo 
we are left with the calculation of 

{C (t)) = A p { -^^- £ J2 e"[t + A ^] . (35) 

P L=2 n=l 

Since Xl(ti) = 2J(mr/L) 2 + 0(1/L A ) it is clear that the dominant contributions to 
fl35|) when t — > oo, are basically contained in large L segments. Introducing the 
(dimensionless) scaling variables 

r = P' s = |> ( 36 ) 

it is straightforward to verify that within the scaling regime £ ^> 1 with Ji ^> 1 
Eq. (|35|) can be written as 

(C( r )) ~ ^ f e~ s £ e- 2w2nV / s2 rfs. (37) 

S ^ n=l 

For rCl, the sum over n can be replaced by an integral, and a result (C(t)) oc t^ 1 / 2 
(of pure-system form) follows. This is the asymptotic behavior on one side of a 
crossover occurring at r ~ 1 . The other side has the more interesting (disorder - 
dominated) asymptotic behavior, obtained by considering t>1. Then the sum is 
dominated by the n — 1 term, so 

(C(r))~f J e'^^ds. (38) 



o 



To obtain the asymptotic behavior of the integral we use a saddle-point expansion 
around the minimum of the exponent, at Sq = (47r 2 r) 1//3 , which becomes exact in 
the scaling limit r ^> 1. After elementary manipulations, this finally yields autocor- 
relations characterized by a universal scaling function U{j) exhibiting a stretched 
exponential decay together with a subdominant power law prefactor. More explic- 
itly, for t > 1 , £ > 1 , 

(C(£,r)> = r'WCr), 



U(r) = aJ^tV* expf-^r^) , (39) 
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where a = (2tt) 2 ^ 3 . This is the central result of this sub-Section, and it corresponds 
to localization when the average cluster size is smaller than the pure diffusion length. 

To check the reliability of this form we address the reader's attention to Fig. 1 
in which autocorrelations are computed directly from Eq. ([35]) . The data collapse 
obtained for different dilutions clearly confirm Eq. fl39|) and provides further evidence 
for the existence of an asymptotic scaling regime characterized by the variable r = 
Jt/l; 2 and showing the crossover at r ~ 1 between pure system diffusive behavior 
(r < 1 ) and localization (r > 1 ). Notice that so long as p ^ 1 , no matter how small 
the dilution is, it induces at sufficiently large time a different dynamics from that 
developed by the regular, purely diffusive case [ Eq. ( |28| ) ] . This contrasting behavior 
comes from the interplay of a pure diffusion length (oc t 1 / 2 ) and the percolation 
length £ (roughly the average segment length) which diverges at the threshold p = 1 . 

B. Percolation cluster 



In higher dimensions the dynamic behavior is much richer by virtue of i) the 
more interesting fractal structure of the percolation geometry at the threshold p c , 
and ii) the existence of the infinite (spanning) cluster above p c . We therefore 
have to distinguish finite cluster contributions to (C(t)) from the infinite cluster 
contribution when p > p c . 

We first consider how such features might affect the Goldstone argument [211]. 
The spin symmetry properties remain as for the pure system. The full rotational 
symmetry by an angle a around an arbitrary spin direction n, evidently leaves the 
evolution operator (Q) invariant, i.e. H = e ~ iaS - n }j e ,<,s '™ . However, the steady 
(ground) states of Eq. @ do not preserve such symmetry ^ z. Because of 
this spontaneous symmetry breaking, in the thermodynamic limit iV — > 00 a low- 
lying band of gapless or Goldstone modes can be expected in the pure system 



and, because of its unlimited size, also on the infinite cluster {p > p c ) irrespective 
of the hopping disorder in H. Moreover, the Goldstone modes are the | A ) states 
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already introduced in Eq. (|19"|) and those with large characteristic scale are ultimately 
responsible for the emergence of a slow asymptotic kinetics [Eqs. (]26f) ,(|ffi) ] . 

Next consider the finite clusters (for p ^ p c ). They are distributed with size 
distribution having characteristic scale £ but nevertheless extending with exponential 
tail to infinite size (see also Section 4). Without these tails the finite size cutoff would 
introduce a gap in the spectrum of the (Goldstone) modes on the ensemble of finite 
clusters. 

The contributions, to dynamic properties such as (C(£, t)) , from excitations on 
finite clusters and on the infinite cluster (for the case p > p c ) can be obtained from 
the density of single levels u(X) or density of states (DOS) involved in Eq. (|25*|). 
Therefore, the averaged autocorrelations can be written as 

oo 

(C(t)) = A p J u(X) e- xt dX . (40) 
o 

The specific form of the DOS for A — > from modes on the finite clusters is 
closely related to the role played by Lifshitz tails ~ e~ x r (r > ) in the energy 



distribution, a typical issue occurring in the presence of disorder |21]| . Consider, for 
example, a square lattice of which a fraction p of bonds are accessible. Below the 
percolation threshold the low energy DOS is characterized by Lifshitz tails which 
are contributed by large regions of connected sites. In the one dimensional case 
(Section 3 A) strings of length L occur with probability p L and contribute low-lying 
ferromagnetic modes with A oc L~ 2 , so the DOS form is cu(X) oc p x 1/2 . In fact, such 
a tail in the DOS is dominant in determining the relaxation dynamics referred to in 
Section 3 A. 

Above the percolation threshold Lifshitz tail modes are also present, but are 
swamped by the power law tail uo(X) ~ A^ 2 " 1 due to low energy excitations on the 
infinite cluster. Thus for p > p c we expect the usual diffusive decay (C(t)) oc t~ d l 2 
at long times. In passing, it is instructive to check the case p = 1. It is well known 
that for an hypercubic lattice the DOS of large wavelength excitations behaves as 



16 



c<j(A) ~ 2{2n) d A^ 2-1 j where Vld is the surface of the ci-dimensional unit sphere. Thus 
after integrating (|40|) , we recover precisely both the asymptotic diffusive tails and 
amplitudes already obtained in Eq. (p8|). 

At p = p c , the Lifshitz tail is again present in uj(X) and swamped by a power 
tail which in this case has a non-trivial exponent d s — 1 where d s is, by defini- 
tion, the spectral dimension and is the ratio of the fractal dimension df and the 
dynamic exponent z (which are respectively the length scaling dimensions for mass 
and frequency): d s = df/z. These exponents have been investigated for the linear 
spin wave problem at the percolation threshold, by direct numerical calculation 
and via scaling relations flID|, [14|j df — d — /3/v, z = 2 + to the exponents for 
percolation conductance (t), correlation length (z/) and density (f3). For all d > 1 , 
d s is very close to (though not exactly equal to) the conjectured Alexander-Orbach 
[23|j value 2/3 (a factor of 2 arises here because the basic equation for ferromagnetic 
spin waves involves frequency squared). 

Hence from Eq. ( f40|) we immediately obtain for the exact finite-density autocor- 
relation function of hard core diffusing particles at large times at the percolation 
threshold (p = p c ) 

(C(t)> oc r d ° , (41) 

where d s = [dv — (3)/ (2u + t — (3) is very close to 2/3 for all d > 1 . 

This is in fact the (t/£ z ) — > scaling limit of the following general form implied 
by scaling considerations applicable when both t and £ are large 

(C(U))*t- d °f(t/¥). (42) 

This form includes the crossover at t ~ £ z , and leads us to conjecture that at p c 
within the percolation cluster, the diffusive behavior is characterized by an anoma- 
lous root mean square displacement R oc , z ~ 3df/2 where df is the fractal 
dimension d — fiv . Eqs. (f4lD and ( ^2|) are the main results of this subsection. 
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4 Conclusions 



The mapping to quantum spin systems exploited in this paper gives a powerful 
approach to the stochastic dynamics of hard core particle systems. Using it we 
have been able to obtain in Section 2 exact dynamical properties (including the 
space and time dependent pair correlation function) for the interacting system at 
any density in terms of corresponding properties for linearized spin wave dynamics. 
The procedure applies equally well to disordered systems, and in addition Goldstone 
arguments still apply there. In the applications made here to disordered systems 
we have chosen to emphasize the diluted case, because of the dramatic effects seen 
there by virtue of the underlying percolation geometry. 

A somewhat surprising result is that hard core particle diffusion on diluted chains 
is quite interesting, any dilution disconnects the chain, so the static diffusion con- 
stant strictly vanishes. Yet finite-time dynamics is non trivial, exhibiting scaling 
behavior and stretched exponential autocorrelation functions in the time domain 
(Section 3 A). The scaling behavior is analogous to that calculated earlier for the ap- 
parently much simpler (but equivalent) problem of ferromagnetic spin wave dynam- 
ics of diluted chains ||24j| . There the calculations were carried out in the frequency - 
wave vector domain, for comparison with inelastic neutron scattering measurements, 
so the stretched exponential time behavior was not identified, but it should also be 
present there. 

For interacting particle dynamics on higher dimensional diluted lattices, by virtue 
of the mapping we have been able to draw on scaling procedures and results for spin 
wave dynamics of diluted systems. This is particularly direct for the autocorrelation 
function of the diffusing hard core particle system, which is given completely in terms 
of a DOS which is the same as for spin wave dynamics. 

Below the percolation threshold the diffusion is entirely on finite clusters. The 
low frequency density of states determining the long time behavior comes from 
large clusters whose occurrence probability is exponentially small in cluster size, 



18 



very much as in the one dimensional case. Above the percolation threshold the 
finite cluster contribution is similar, but it is now swamped by the infinite (span- 
ning) cluster contribution which, well away from the threshold gives a behavior like 
that for the pure system. At and near the threshold (p ~ p c ) the diverging scale 
of the percolation geometry induces scaling behavior for the long time dynamics. 
Here the DOS has an anomalous power law dependence on frequency A; the asso- 
ciated exponent has been obtained via the relationship of the DOS to that for the 
linear spin wave system and further relationships to percolation processes including 
(static linear) conductance; the exponent in the DOS power law is given in terms of 
t, p, u, the exponents dependence on \p — p c \ of percolation conductance, density 
and correlation length. The consequent exponent in the power law time decay of the 
autocorrelation function at p c turns out to be close to 2/3 for any d > 1 (in accord 
with the approximate Alexander-Orbach conjecture ). The scaling behavior of 
the space and time dependent correlation function depends on the DOS exponent 
and also on the length scaling exponent for the frequency A . This so called dy- 
namic critical exponent is also related to t, (5 v , and hence we can give the power 
law dependence of the characteristic diffusion length on time. This generalization 
also allows discussion of the crossover due to the competition between characteristic 
diffusion length and percolation correlation length near the transition. 

Clearly, further extensions are desirable. For the dilute systems the most obvious 
ones are the calculation of space and time dependent correlation and crossover func- 
tions, and also probability distributions rather than just averages. In addition sim- 
ulation confirmation of the analytic predictions or experimental comparisons would 
be desirable. As regards more general disorder, other types of substitutional disorder 
could clearly be treated by the methods we have used. A more challenging extension 
is to the biased diffusion of hard core particles: much has been done exploiting the 
quantum spin mapping for the biased case along with adsorption-desorption pro- 
cesses [f| . But biased hard core particle diffusion on disordered lattices is at present 
beyond the exact techniques exploited here. 
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Figure 1: Stretched exponential decay and scaling regime for autocorrelation func- 
tions of diluted chains, r (here denoted as T), is the time t divided by the square 
of the percolation correlation length £ . Following Eq. (38) in the text, F(t) is here 

taken as (C(£, t 1 / 6 \J^f-\. Solid lines denote different degrees of dilution 

p = e~ 1 ^ . The uppermost curve at the left corresponds to p — 0.6, and the lower 
curves to p — 0.5, 0.4, 0.3, 0.2 and 0.1, in descending order. The slope of the dashed 
line is -|(27r) 2/3 log 10 e. 
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